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The atmosphere of Jupiter has east-west zonal jets that alternate as a 
function of latitude as tracked by cloud motions at tropospheric levels. 
Above and below the cold tropopause at -100 mbar, the equatorial 
atmosphere is covered by hazes at levels where thermal infrared 
observations used to characterize the dynamics of the stratosphere lose 
part of their sensitivity. James Webb Space Telescope observations of Jupiter 
inJuly 2022 show these hazes in higher detail than ever before and reveal the 
presence of an intense (140 m s”) equatorial jet at 100-200 mbar (70 ms? 
faster than the zonal winds at the cloud level) that is confined to +3° of the 
equator and is located below stratospheric thermal oscillations that extend 
at least from 0.1 to 40 mbar and repeat in multiyear cycles. This suggests that 
the new jet is a deep part of Jupiter’s Equatorial Stratospheric Oscillation and 
may therefore vary in strength over time. 


12,13 


Jupiter and Saturn atmospheres have zonal jets with broad prograde relationship”, reveal the presence of intense stratospheric eastward 


winds in their equatorial regions’”. The winds are best observed at 
cloud-top levels but extend 3,000 km deeper inJupiter** and 8,000 km 
inSaturn>°. Above the cloud level, wind speeds decay with increasing alti- 
tude tonearly zero values at 20 mbar in most of Jupiter’s mid-latitudes’®. 
However, this pattern of vertical decay is spectacularly broken in Saturn’s 
equator, where an intense (400 m s”) and narrow (5° width) jet is observed 
tracking the motions of high-altitude hazes at the 50-60 mbar level”, 

In the equatorial regions of both planets, temperature retrievals 
through infrared spectroscopy, in combination with the thermal wind 


and westward winds alternating in direction as a function of altitude 
in the 0.1-40.0 mbar region®"*. The perturbations to the equatorial 
stratospheric temperatures and inferred thermal winds vary periodi- 
cally on multiyear timescales’, sharing characteristics with similar 
phenomena on Earth, suchas the quasi-biennial oscillation (QBO) and 
semi-annual oscillation (SAO)”””. Like their terrestrial counterparts, 
Jupiter and Saturn’s equatorial stratospheric oscillations form a ver- 
tical pattern of temperature and wind perturbations that gradually 
descends over time’? ”. Jupiter’s equatorial stratospheric oscillation 
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Fig. 1| NIRCam images of Jupiter. Images are shown after combining individual 
frames in the same filter which differed in their acquisition time by a few 
minutes. These frames were first corrected for the rotation of the planet. 

A high-pass filter was applied to enhance the contrast. A linear-stretch version of 
this figure without high-pass processing is available as Supplementary Fig. 1. 
Images in the F164N filter were obtained in subarrays and mosaicked into 
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a single image compensating for the rotational effects of the 16 individual 
frames that were used. Images in the F164N and F212N filters have a pixel scale 
of 0.03087 arcsec pixel (roughly 99 km pixel” on Jupiter’s disc). Images in the 
F335M, F360M and F405N filters have a pixel scale of 0.06303 arcsec pixel * 
(203 km pixel”). 


(JESO) has a variable period of 3.9 to 5.7 years with the variability linked 
to the occasional development of large-scale convective perturba- 
tions”. Saturn’s equatorial stratospheric temperatures oscillate 
quasi-periodically in a 15-year cycle (that is, half a Saturn year)", but 
this period can also be disrupted by large convective storms develop- 
ing at other latitudes”. Variations of temperatures in Jupiter’s upper 
troposphere have been observed down to 300 mbar witha periodicity 
at the equator of 8 years and an anticorrelated relation with the JESO at 
near-equatorial latitudes”. The complex phenomenology at the equa- 
torial regions of Jupiter and Saturn points to connections between the 
stratospheric and tropospheric dynamics that are not well understood. 

Elevated hazes cover Jupiter’s equatorial region around and above 
the 200 mbar level” ”. These hazes are variable in opacity and albedo, 
and their brightness in the methane absorption band at 890 nm can 
vary upto 30% over different years”. The equatorial hazes are situated 
well above the tropospheric clouds, where winds are well-characterized 
by trackable cloud features, and well below the stratospheric region, 
where temperature retrievals from thermal infrared observations allow 
deducing winds through the thermal wind analysis. However, the hazes 
appear mostly featureless, preventing high-quality measurements of 
winds at their altitude”’*’. Ultraviolet (UV) images acquired by the Cas- 
sini spacecraft, representative of the 500 mbar level, showed that right 
at the equator about 4-5 km above the cloud tops at around 600 mbar, 
zonal winds increase slightly (10 m s”) with increasing altitude”, but 


those measurements had large uncertainties due to the low contrast 
in UV images. 


Observations 

We observed Jupiter’s atmosphere on 27 July 2022 using the Near 
Infrared Camera (NIRCam) instrument** on the James Webb Space 
Telescope (JWST) as part of the Jupiter System Early Release Science 
Program no. 1373 (ref. 35). Because Jupiter is a bright target for JWST, 
we selected filters centred in strong atmospheric methane absorp- 
tion bands (F164N, F335M, F360M and F405N) and in the strong H,-H, 
and H,-Hecollision-induced absorption band at 2.12 um (F212N), and 
obtained full-disk images of the planet. We repeated our observations 
after one planetary rotation (about 10 h) to observe the same hemi- 
sphere twice and examine temporal changes and dynamics, with the 
exception of filters F360M and F405N, which were only used once. 
Figure 1 shows NIRCam images processed with a high-pass spatial filter 
to emphasize fine-scale details (Supplementary Table 1 gives details 
about the images and Supplementary Fig. 1 shows these images without 
high-pass filtering). The Methods provides details about the image 
acquisition and processing techniques. 

Ground-based and Hubble Space Telescope (HST) observations 
of Jupiter at these wavelengths generally result in images with very 
low contrast in which most fine-scale atmospheric details cannot be 
resolved*’. Our JWST observations clearly resolve fine-scale details all 
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Fig. 2 | Pressure at which two-way vertical optical depth t reaches 1 (black 
line) for a model atmosphere with no aerosols. The model atmosphere with no 
aerosols is based on ref. 37. The yellow shaded region shows the two-way vertical 
optical depths from 1 to 2. The grey-blue shaded area shows the two-way vertical 
optical depths from 2 to 5. The curves are compared with system throughputs 
for different filters. Filters marked with * are those for which winds have been 
obtained in the past (F631N, and the CB2 filter on the Cassini/ISS), or in this work 


(F164N, F212N, F335M). For each filter, the mean pressure level sensed (T = 1) is 
given at the top, and the range of levels that can be sensed with optical depths of 
2-5 is given next to each filter in the figure. A is the light wavelength measured in 
pum. a, Selection of filters in which atmospheric dynamics have been studied in 
the past (F631N filter in HST/WFC3 instrument, the CB2 filter on Cassini/ISS and 
the FQ889N filter in HST/WFC3 centred at a strong methane absorption band. 
b, Filters used in our JWST/NIRCam observations. 


over the planet at all wavelengths. Thus, the image quality provided by 
JWST enables new capabilities to study Jupiter’s atmospheric dynamics 
at the altitudes sampled. 

The JWST filters in which we observed Jupiter are sensitive to a 
mixture of tropospheric and stratospheric hazes (Fig. 1; filters F335M, 
F360M and F405N are also sensitive to H}* emissions from polar auro- 
ras which will be discussed in a separate paper). Figure 2 shows the 
sensitivity of JWST/NIRCam filters to different altitudes for a simpli- 
fied model atmosphere that assumes no aerosols and considers only 
gaseous absorptions and Rayleigh scattering”. We show the two-way 
penetration depth of solar radiation in Jupiter’s atmosphere consider- 
ing optical depths from t= 1to T= 5. Sensing levels by these filters will 
depend on the structure of clouds and hazes and will be different in 
different locations of the planet. For each filter, regions with clouds 
or hazes located at altitudes with low optical depths will appear bright 
(an example is the Great Red Spot (GRS)) (Fig. 1). In filters with strong 
atmospheric absorption bands, regions with tropospheric clouds 
without substantial haze above them will appear dark (like the area 
around the GRS and most locations over the planet in Fig. 1). 

Inthe filters F164N, F212N, F335M and F360M, unenhanced images 
in Supplementary Fig. 1 show large brightness differences between 
the equatorial zone (EZ), the GRS, the polar areas and the rest of the 
planet. Inthe F405N filter, which is not affected by notable absorption 
from gases, the contributions of low optical depth small particle hazes 
are greatly diminished, and the tropospheric levels where the clouds 
become optically thick are sensed. Comparing these multiwavelength 
images, we see that outside of the EZ, the GRS and the polar areas, all 
JWST images show the same features sensing the main cloud deck in 
the troposphere (500-700 mbar). 

The EZ inthe F164N, F212N, F335M and F360M images is roughly 5, 
9,4. and 2times brighter respectively than the mid-latitudes. The high 
brightness difference of the EZ from the rest of the atmosphere in the 
F212N and F164N filters suggests that in both filters we are observing 
the low-opacity hazes at elevated altitudes close to T = 1, while in filters 
sensitive both to the lower stratosphere and the upper troposphere 
(F335M and F360M), the contrast between elevated and deep features 
is reduced. At the EZ, these two filters have contributions from haze 
layers above those sensed with the F212N and F164N filters mixed with 
deeper features in the troposphere. The approximate levels sensed in 


the EZ in our observations are summarized in Table 1, where they are 
compared with levels sensed in the past to investigate cloud motions. 

JWST/NIRCam images show an amount of detail in the equatorial 
hazes (Fig. 3) that largely surpasses those of any previous images of 
Jupiter at these wavelengths. The 2.12 um images show crisp features 
(highlighted in Fig. 3) that have bright counterparts in all other JWST 
filters and seem to be associated with bright and compact clouds inthe 
lower troposphere as observed in the F405N image. These features may 
represent either convective activity, or a vertical connection between 
tropospheric clouds and features located at least one scale height 
above the main cloud deck. 

Supplementary Fig. 2 shows HST images and maps of the equato- 
rial region obtained on 28 July 2022 two Jupiter rotations after the )WST 
images shown in Fig. 3. Some of the bright cloud systems in the F405N 
images show convective activity in these HST images of the cloud level. 
Time-resolved ground-based observations of Jupiter obtained through 
July 2022 and available at the PVOL database (http://pvol2.ehu.eus) and 
ALPO Japan website (https://alpo-j.sakura.ne.jp/indexE.htm) confirm 
that the planet’s atmosphere was relatively quiescent, with no major 
convective outbreaks. 


Wind results 
Wetracked the motions of atmospheric features observed in the F164N, 
F212N and F335M images by comparing map-projected images of the 
planet separated in time by one planetary rotation. We used image 
correlation software designed for cloud-tracking wind analysis** *° 
and excluded from the analysis regions covered by large atmospheric 
systems like the GRS (Methods). Figure 4 shows our wind results. 
Outside of the EZ, our new JWST-based zonal wind profiles are 
very similar to prior measurements of Jupiter’s zonal winds at the cloud 
tops*”*"? (Fig. 4a). This is because out of the EZ, the cloud systems 
tracked in NIRCam images are the same as those visible in HST and 
ground-based images (Supplementary Fig. 2), that is, the absence of 
elevated aerosols across most of the planet allows these filters to sense 
the lower troposphere. The quality of the wind retrievals at high lati- 
tudes (Fig. 4a) is progressively deteriorated by the increasing opacity 
of elevated subpolar hazes that obstruct views of the tropospheric 
clouds. Throughout the EZ, at latitudes +10° (all latitudes in this paper 
are planetographic), the derived zonal winds represent the motions 
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Table 1| Approximate sensing levels in the EZ when taking 
into account clouds and hazes 


Source Filter A. (um) Pressure Primary sensed aerosol 
level layers EZ 
(mbar) 
UV1 0.264 315-500* Equatorial haze 
HST/WFC3 and F631N 0.630 500-600" Tropospheric clouds 
Cassini/ISS CB2 0750 500-600" Tropospheric clouds 
FQ889N 0.889 200-400" Equatorial haze 
F212N 2120 ~240 Equatorial haze 
FI64N 1.644 ~200 Equatorial haze 
F335M 3.365 50-500 Equatorial haze+ 
JWST/NIRCam tropospheric clouds 
F360M 3.621 100-500 Equatorial 
haze+tropospheric 
clouds 
F405N 4.055 500-600 Tropospheric clouds 


Àe is the pivot wavelength for each filter”. CB2, Continuum Band 2; ISS, Imaging Science 
Subsystem; WFC3, Wide Field Camera 3. * From ref. 33. ° From ref. 29. 


of the upper hazes, which show a very different behaviour from the 
motions at the main cloud layer. A distinct central narrow and intense 
jet is observed at the three wavelengths where JWST obtained data on 
two consecutive rotations of the planet (F164N, F212N and F335M). 
Outside of the central narrow jet, but still inside the EZ, JWST images 
show a decrease in the zonal winds at higher altitudes with respect to the 
winds at the cloud tops. This decrease in the wind intensity with altitude 
has some minor differences between the northern edge of the EZ and 
the southern edge. The latitudinal asymmetry in the winds may relate to 
different haze opacities north and south of the equator.JWST observed 
higher haze opacity in the southern branch of the EZ, as illustrated in 
Fig. 3, but we note that previous studies have shown that the latitudinal 
asymmetry in haze distribution is variable on timescales of a few years”. 

In the F212N observations, zonal winds in the EZ from the cloud 
correlation analysis have lower dispersions than in the F164N and 
F335M filters. In the F212N images, details are better contrasted and 
individual cloud features can be visually identified and tracked after a 
planetary rotation (Supplementary Fig. 3 displays a selection of those 
features, blue dots in Fig. 4b depict the velocity of the full sample of 
visually identified cloud-tracking details). However, many of the crisp 
features highlighted in Fig. 3 have no counterparts one rotation later, or 
appear very distorted on the next planetary rotation, which is consist- 
ent with the presence of strong wind shears. 

The wind measurements in the EZ obtained in the F164N and F335M 
filters have larger dispersions than those fromthe F212N images (Fig. 4c). 
This is caused by a combination of different factors. First, the wind 
measurements in F164N images were obtained in subarrays with nar- 
rower views of the EZ. This results in smaller sizes of the regions where 
cloud correlation can be calculated. Secondly, the F335M images have 
halfthe spatial resolution of the F164N and F212N images and are also 
subject to higher noise levels. Finally, the F164N and F335M images are 
both sensitive to a wider vertical column of the atmosphere than the 
F212N images. Thus, we cannot discount the possibility that part of the 
wind dispersion in these two filters is created by tracking features at 
different altitudes through the vertical range of levels sensed by these 
filters and affected by vertical wind shear. 

Our measurements show that the zonal winds at the equator 
increase in height from the cloud tops at around 600 mbar inthe EZ”, 
where the zonal wind profile has a local minimum of 70 ms“, to the 
tropopause at 100-200 mbar, where the zonal wind profile has alocal 
maximum of 140 ms inthe F335M filter. The different wind values at 
the equator obtained in F212N (105 +16 m s”), F164N (115 + 35 m s”) and 


F335M (140 + 20 m s”) images suggest the presence of a strong vertical 
wind shear within the haze levels, with the most intense winds observed 
inthe F335M images most likely sensing the highest levels (Supplemen- 
tary Fig. 4). Furthermore, UV-measurements of zonal winds during the 
Cassini flyby in 2000 found zonal winds of 82 + 11 m s” at 500 mbar 
(ref. 33), supporting the hypothesis that most of the vertical wind shear 
in this central equatorial jet is concentrated through the hazes. Right 
at the equator, we estimate mean vertical wind shears from the cloud 
tops to the hazes of 40-70 ms" per scale height (-20 km), with more 
uncertain wind shears of 10-110 m s” per scale height within the hazes 
(Supplementary Fig. 4). Outside the central jet at +3°, but still in the EZ 
(+10°), the winds decrease in height, with the strongest shear at +7°, 
where the zonal winds decrease from the cloud tops to the equatorial 
hazes by 45-60 ms”. This results in a mean vertical wind shear of -25 
to-60 ms ‘per scale height, which is comparable but opposite in sign 
to the vertical shear at the equator. 


Discussion 
The new equatorial jet near Jupiter’s tropopause makes Jupiter’s and 
Saturn’s equatorial circulations remarkably similar (Fig. 5). Both planets 
harbour a doubly peaked equatorial jet at the main cloud level observ- 
able in visible wavelengths. Our new JWST observations show that 
Jupiter has in addition a sharply peaked jet right at the equator (+3°) 
inthe hazes above the main cloud deck. The wind speed in the narrow 
jet increases with altitude up to the tropopause, while it decays with 
altitude outside of the central jet. A similar wind structure was found 
previously on Saturn’s equator, except that the central jet inJupiter 
increases in altitude much more strongly than on Saturn. Jupiter’s cen- 
tral equatorial jet is located in a region of the EZ that appears brighter 
than the rest of the EZ at 889 nm and at F335M (Fig. 5). On Jupiter and 
Saturn, the zonal wind profiles at cloud level have two local maxima near 
+10° inJupiter and +15° in Saturn that move poleward with increasing 
height. Besides these similarities, there are also differences between 
the equatorial winds in both planets. Jupiter’s EZ is narrower in latitude 
than Saturn’s EZ and the observed meteorology at cloud level is very 
different”. In addition, in Jupiter the increase in wind intensity with 
altitude at the equator (40-70 m s” per scale height; see also Supple- 
mentary Fig. 4) is notably larger than the one in Saturn (8 ms“ per scale 
height in the Cassini epoch)". However, these vertical shears could be 
variable in time if the jets at the upper hazes followtemporal variations 
similar to the oscillations in stratospheric temperatures and winds. 
The equatorial stratospheric thermal oscillations and asso- 
ciated winds in Jupiter and Saturn also share similarities and 
differences? 9-5, Among the latter, the different periods of JESO 
(3.9-5.7 years”) and Saturn’s SAO (15 years”) argue for different 
dynamical origins. The half Saturn year cycle of Saturn’s SAO suggests 
a seasonal effect”, while the time periods and phenomenology of the 
JESO cannot be related with seasons. JESO is reminiscent of the Earth’s 
QBO”®*, a downward propagating easterly and westerly wind pattern 
that affects the 5-50 mbar portion of the Earth’s atmosphere. Although 
the QBO phenomenology on Earth is not fully understood, gravity 
waves created by convection in the lower troposphere are known to 
playa major role releasing energy and momentum inthe Earth’s strato- 
sphere”. Similar mechanisms, suchas the forcing created by gravity 
waves originating from convection in the troposphere, have been 
proposed for Jupiter’s thermal oscillations“. Our JWST observation in 
F405N in Figs. 1 and 3, and the context views of the lower atmosphere in 
Supplementary Fig. 2, show bright cloud systems of possible convective 
nature with small-scale bright systems observable in the upper hazes. 
OurJWST observations provide new insights into the troposphere- 
stratosphere interactions in Jupiter, and can inform us about the depth 
of penetration into the troposphere of the equatorial oscillation. On 
Earth, the QBO does not penetrate below the tropopause, but on 
Jupiter, the unexpected presence of this intense and elevated equa- 
torial jet, which is not reproduced in existing simulations, suggests 
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Fig. 3 | Multiwavelength maps of the equatorial region. a-e, Images in a-d upper hazes. Purple arrows indicate bright details observable in the F405N image 
were obtained close in time and the image in e was obtained after one planetary with counterparts only in some of the images of the upper hazes. Yellow arrows 
rotation (Supplementary Table 1), using filters F164N (a), F212N (b), F335M show a convective storm also active in HST observations taken one rotation 

(c), F360M (d) and F405N (e). All images have been processed with a high-pass later (Supplementary Fig. 2). The green arrow in e shows the location of a deep 
spatial filter to show the smallest-scale details. Boxes show bright features in cloud system that developed into a convective storm one Jupiter rotation later as 
the lower troposphere observed in the F405N image with counterparts in the observed in HST images. 
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Fig. 4 | Jupiter’s zonal wind profile. Zonal winds from JWST images are 
compared with zonal winds measured on Cassini UV images”, on continuum 
filters at cloud level“, and on HST images acquired in 2019 and analysed following 
the procedure in ref. 39. a, Mean zonal winds. b, Zonal winds at the equator. 

A black line shows the average of Cassini data from the year 2000 and HST data 
from 2019. The blue line is the F212N zonal mean winds. Grey symbols correspond 
to measurements of F212N images using the cloud correlation software. Blue 
circles correspond to cloud features visually tracked in F212N images. Mean error 
bars from twice the standard deviation (20) are shown as segments for each wind 
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profile in a and for the +10° latitudes in b. c, Zonal winds in the F164N, F212N and 
F335M filters in the equator compared with cloud-top winds and winds measured 
on Cassini UV images. Dots are individual measurements from the correlation 
software and lines are mean zonal winds. The inset shows the standard deviation 
of zonal winds in each filter considering latitudinal bins of width 0.3°, except for 
the Cassini UV data, which used latitudinal bins of 1°. d, Wind difference between 
the different filters sampling the hazes and the clouds tops. A dotted grey line 
shows the zero value. 


the stratospheric equatorial oscillations penetrate into the upper 
troposphere in an unanticipated way. The stratospheric oscillations 
are variable in both Jupiter and Saturn, and the time variability of the 
jet at the base of the stratosphere remains uncertain given that we have 
only observed Jupiter’s lower stratospheric jet during a single epoch. 
Thermal maps of Jupiter acquired in mid-2019 (ref. 20), extrapolated in 
time to 2022 using the descent rate of the JESO observed during its last 
cycle, suggest that the -40 mbar level had a warm anomaly during our 
JWST observations. This would be consistent with winds decreasing at 
greater depths as observed in our JWST images. This extrapolation in 


time also suggests a change in the sign of the thermal perturbation at 
40 mbar in 2023-2024, allowing an examination of the time-variable 
component of the equatorial jet. 

A combined view of Jupiter’s equatorial dynamics from the cloud 
tops to the upper stratosphere can be obtained by combining observa- 
tions in the visible (500-700 mbar), near-infrared (100-200 mbar), 
thermal infrared (1-40 mbar), and Doppler winds from millimeter data 
(-4 mbar), where intense equatorial zonal winds have been observed**”’. 
At 500-700 mbar, wind measurements at visible wavelengths hint at 
variability (<10 ms“) with possible periods of 5-7 years*®, while at 
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Fig. 5 | Sketch of Jupiter and Saturn equatorial dynamics. a, Jupiter’s EZ and 
equatorial winds from Cassini (green line cloud tops, white dots UV), HST (blue) 
and JWST F335M wind profile (red). Background images show the morphology 
ofthe EZ in different wavelengths. b, Saturn’s EZ and equatorial winds from 
Voyager (green line), Cassini deeply penetrating filters (blue) and Cassini filters 
sampling the 889-nm methane absorption band (red). Background images show 
Saturn’s EZ in different wavelengths observed in 2022 by HST, with the grey 


T(K) 


T(K) 


region covered by the rings. c, Sketch of the different cloud layers and levels 
where zonal winds are observed (circles). Both the cloud-top winds and the winds 
at the hazes are eastward, but winds in the upper stratosphere alternate eastward 
or westward depending on the altitude level and the phase of the thermal 
oscillation. P, atmospheric pressure measured in millibars; T, atmospheric 
temperature measured in Kelvin. 


stratospheric levels the periodicity in the temperature oscillations and 
winds is 4-6 years””™. In the near-infrared windows in which the upper 
hazes are observable at intermediate levels between the cloud tops and 
the stratosphere, only new observations by JWST can determine the 
atmospheric circulation and its expected variability. The equatorial 
dynamics of Jupiter and Saturn are complex, with similarities and dif- 
ferences that can guide us to unveil the various mechanisms that shape 
atmosphere dynamics at low latitudes in fast-rotating giant planets. 


Methods 

Image acquisition 

Jupiter represents a bright target for JWST. We explored the various 
filters and strategies in which Jupiter could be observed with NIRCam 
without saturating the detectors”. Image saturation was avoided by 
exposing the planet in filters centred in strong absorption bands of the 
Jovian atmosphere (F212N, F335M) and using the shortest integration 
times possible. In filters where the planet is brighter (F164N, F360M 
and F405N) we used subarrays (SUB640) that reduce the minimum 
exposure time. Inall cases we used rapid readout patterns. We designed 
a pointing sequence in which gaps between the four NIRCam B detec- 
tors can be recovered with an appropriate dither pattern using four 
exposures separated in time by only a few minutes. Images obtained 
using a full array configuration were acquired with a dither pattern of 
three positions. 


Image processing and composition 

We used the NIRCam calibration processing pipeline” to produce cali- 
brated and geometrically corrected versions of the images with differ- 
ent exposure levels from the original uncalibrated images. NIRCam 
detectors were read out non-destructively three times for each image, 
conserving the charge in each pixel and forming different versions of 
the same image (groups) with different accumulated exposures. We 
used a combination of one, two and three groups in images obtained 
in full arrays, and versions of the images from only two and three 
groups in images obtained using subarrays to remove saturated areas, 
produce sharp images without rotational smearing, and obtain high 
signal-to-noise ratio through our images. We used an adaptive median 
filter to remove bad pixels identified by the pipeline. We navigated each 
individual image and subframe separately comparing the position of 
the planet witha synthetic longitude-latitude grid using the WinJupos 
software due to its flexible use. These navigations were compared with 
navigation obtained with python code using World Coordinate System 
information from the telescope pointing, finding good agreement. 
Images obtained with small time separations of a few minutes were com- 
bined with WinJupos correcting the rotation of the planet and produc- 
ing higher signal-to-noise-ratio images. These were high-pass filtered 
to enhance the contrast of small-scale faint features in the planet and 
reduce limb-darkening effects. The images were map-projected, over- 
sampling the initial resolution of the images, that is the resolution used 
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inthe maps was greater than that of the original images, limb-darkened 
corrected and further enhanced with high-pass filters. 


Winds analysis 

Maps of the planet were compared using the PICV3 software™?. This 
software performs one- or two-dimensional correlations over two 
images using boxes whose size can be configured and adapted to the 
contrast of the region of the image where the analysis is run. Winds 
were measured over the final composed images except in the F164N 
images, where image pairs over the same B detector were compared 
individually producing 16 individual data sets that were later merged 
into a single zonal wind profile. PICV3 has a graphical user interface that 
allows selection of specific regions of the images. We avoided regions 
with large-scale features and examined the equatorial region with long 
boxes with typical longitudinal sizes of 20° and latitudinal sizes of 1°. The 
F164N images were analysed using 16 image pairs from the combination 
of 4 subarrays and 3 dither positions combining the results from the 
individual zonal wind profiles ina single zonal wind profile in that filter. 


Data availability 

Source data for this paper are available at: https://github.com/ 
JWSTGiantPlanets/Jupiter-Atmosphere-NIRCAM. This dataset is asso- 
ciated with the following https://doi.org/10.5281/zenodo.8199030. 
JWST and HST data used in this paper were obtained from the Mikulski 
Archive for Space Telescopes at the Space Telescope Science Institute 
(https://archive.stsci.edu/), which is operated by the Association of 
Universities for Research in Astronomy, Inc., under NASA contract 
NAS 5-03127 for the JWST and under NASA contract NAS 5-26555 for the 
HST. JWST ERS 1373 observations of Jupiter’s atmosphere used in this 
publication are available at https://doi.org/10.17909/vc7n-p631. HST 
context observations obtained on 28 July 2022 are available at: https:// 
doi.org/10.17909/rnp9-gj41. HST/OPAL maps of Jupiter from 2019 and 
Saturn from 2022 are available at https://doi.org/10.17909/T9G593. 
Amateur observations used to provide temporal context to our JWST 
observations can be downloaded from the PVOL (http://pvol2.ehu.eus) 
and ALPO Japan (https://alpo-j.sakura.ne.jp/indexE.htm) databases. 
Source data are provided with this paper. 


Code availability 

JWST/NIRCAM calibration software is available at https://github.com/ 
spacetelescope/jwst and is documented online at https://jwst-docs. 
stsci.edu/jwst-science-calibration-pipeline-overview. WinJupos soft- 
ware is available at http://jupos.org/gh/download.htm. The image 
correlation software used in this study (PICV3) is available at https:// 
doi.org/10.5281/zenodo0.4312675. 
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